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Abstract 

A Bayesian non-parametric approach for efficient risk management is pro- 
posed. A dynamic model is considered where optimal portfolio weights and 
hedging ratios are adjusted at each period. The covariance matrix of the 
returns is described using an asymmetric MGARCH model. Restrictive para- 
metric assumptions for the errors are avoided by relying on Bayesian non- 
parametric methods, which allow for a better evaluation of the uncertainty in 
financial decisions. Illustrative risk management problems using real data are 
solved. Significant differences in posterior distributions of the optimal weights 
and ratios are obtained arising from different assumptions for the errors in the 
time series model. 
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1 Introduction 



The first cases of managing risk can be traced as far as to the ancient world, see 



McNeil et al. (2005). Ever since Markowitz (1952) introduced the mean-variance 



approach, the area of financial risk management has advanced immensely. Portfolio 
optimization and hedging are just few out of many risk management problems, how- 
ever, in this volatile world, they are as relevant as ever for today's investor. While 
some individuals seek personal gain, others try to ensure stability and reduce risk. 
Both of these goals can be achieved via efficient allocation and protection of their as- 
sets. For more on modern portfolio theory, hedging and risk management in general 
see |Korn and KoTn] (pOOlj), |Elton et al.\ ([2003|, |McNeil et al] pOOo), |Priei| pOOTl), 



Kwok (2008) and Hull (2012), among others. These papers show that in order to 



determine optimal portfolio weights or hedging ratios, it is very important to find 
appropriate financial models which describe adequately the individual variability of 
the assets and their correlations or mutual dependence. 

Traditionally, risk management problems have assumed a time-invariant rela- 
tionship structure. However, ever since Engle (1982) showed the existence of time- 
varying variances, the standard approach of using constant unconditional correla- 
tions and covariances is being debated in the financial literature. The overwhelming 
empirical evidence shows the advantages of employing a time-varying approach, see 



Rossi and Zucca (2002), Yang and Allen (2004[), Giamouridis and Vrontos (2007), 



Lien (2009), Liu et al. (2010), Lee and Lee 



2012[) and Basak and Chabakauri (2012), 



among others. ARCH-family models, first introduced by Engle (1982) and then gen- 
eralized by Bollerslev (1986), without doubt, are the most researched and used in 



practice to explain time- varying volatilities, see also Bollerslev et al. (1992), Boiler- 



slev et al. 


(1994 


), 


Engle 


(2002b 


), 


Terasvirta 


(2009 


) and 


Tsay 


(2010) 



with multivariate time series, one must also take into consideration the mutual de- 
pendence between returns. Correlations can also exhibit some stylized features, such 



as persistence and asymmetry. For multivariate GARCH (MGARCH), see Bauwens 



et al. (2006), Silvennoinen and Terasvirta (2009) and Tsay (2010). In this paper we 



will consider a very general multivariate GARCH model which accounts for both 
types of asymmetries: in individual volatilities and in conditional correlations. 

Whichever GARCH-type model is chosen, the distribution of the returns de- 
pends on the distributional assumptions for the error term. It is well known, that 
every prediction, in order to be useful, has to come with a certain precision measure- 
ment. In this way the agent can know the uncertainty of the risk she is facing. Distri- 
butional assumptions permit to quantify this uncertainty about the future. However, 
the traditional premises of Normal or Student-t distributions may be rather restric- 
tive. Alternatively, in this paper, we propose a Bayesian non-parametric approach 
avoiding the specification of a particular parametric distribution for the return inno- 
vations. More specifically, we consider a Dirichlet Process Mixture Model (DPM), 
firstly introduced by Antoniak (1974), with a Gaussian base distribution. This is 
a very flexible model that can be viewed as an infinite mixture of Gaussian dis- 
tributions which includes, among others, the Gaussian, Student-t, logistic, double 
exponential, Cauchy and generalized hyperbolic distributions, among others. 

The Bayesian approach also helps to deal with parameter uncertainty in port- 
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folio decision problems, see e.g. Jorion (1986), Greyserman et al. (2006), Avramov 



and Zhou (2010) and Kang (2011), among others. This is in contrast with the usual 



maximum likelihood estimation approach, which assumes a "certainty equivalence" 
viewpoint, where the sample estimates are treated as the true values, which is not 



always correct and has been criticized in a number of papers. As noted by |Jorion 
(1986), this estimation error can gravely distort optimal portfolio selection. In this 
paper, we propose a Bayesian method which provides the posterior distributions of 
the one-step-ahead optimal portfolio weights and hedging ratios, which are more in- 
formative than simple point estimates. In particular, using the proposed approach, 
it is possible to obtain Bayesian credible intervals for the optimal portfolio weights 
and hedging ratios. Also, as seen in Ardia and Hoogerheide (2010), the Bayesian 
inference provides some other advantages over the classical maximum likelihood 
techniques. For example, it is easy to incorporate via priors complicated positivity 
constraints on the parameters to ensure positive variance and covariance station- 
arity. Additionally, it is possible to approximate the posterior distribution of any 
other non-linear function of the parameters, as will be done for the optimal portfolio 
weights and hedging ratios. Moreover, the results are reliable even for finite samples. 
And finally, the models we wish to compare do not necessarily have to be nested. 

Therefore, the main contribution of this work is the application of Bayesian 
non-parametric techniques in portfolio decision problems and exploration of the dif- 
ferences in uncertainty between the proposed approach and conventional restrictive 
distributional assumptions. Our objective is to provide a more realistic evaluation of 
risk of financial decisions. More specifically, in this paper we solve time- varying port- 
folio allocation and hedging problems using the multivariate GARCH specification 
of Cappiello et al. (2006), combined with ideas of Hafner and Franses (2009), and 
univariate GJR-GARCH model of Glosten et al. (1993) for the individual volatili- 
ties. For the errors, we assume a Bayesian non-parametric model based on the class 
of multivariate scale Gaussian mixtures where the scale mixing distribution follows 
a Dirichlet Process (DP) prior (Ferguson, 1973), leading to a DPM model. 

The outline of the paper is as follows: Section [2] introduces the static and time- 
varying portfolio optimization and hedging approaches. Section [3] describes the 
model, inference and prediction from a Bayesian perspective. Section [4] presents a 
short simulation study. Section [5] illustrates the proposed approach using two real 
data examples. Finally, Section [6] concludes. 



2 Portfolio Decisions 

In this section, we first introduce Global Minimum Variance (GMV) portfolio and 
solve portfolio allocation problem by maximizing agent's utility. Then, in order to 
protect the portfolio from market risk, we solve a hedging with futures problem. 
The financial applications are presented using static and time-varying approaches 
and comparing different distributional assumptions in a Bayesian context. 
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2.1 Portfolio Allocation and Hedging: Static Approach 

The main objective of diversification is to increase investor's utility and reduce 



her exposure to risk. See Markowitz (1952) and Merton (1972) for some classical 



portfolio optimization references. In this paper we consider the cases where the 
investor maximizes her expected utility and minimizes the portfolio variance. The 
GMV portfolio can be found at the very peak of the efficient frontier. For the utility, 
assume quadratic preferences. The following are both optimization problems: 

p\j = argmaxE \Uipf, rf)] = argmaxEfrf] — — Var[rf] : p'Ik = 1, 
p p 2 

Pgmv = arg minVar [rf] : p'l K = 1, 
v 

where p is the weight vector, 7 is the risk- aversion coefficient, representing a trade- 
off between the expected return and risk, Ik is a K- vector of ones and rf is a vector 
of portfolio returns. Portfolio is composed of assets, where rf = p'r t and r t is a 
K x 1 vector of asset returns, such that E[rJ = [i and Cov[r t ] = E. The following 
are the closed-form solutions for both portfolios: 



, 1 / ! ^- 1 l K - 1 . \ 

Pu = i^ vjFTk ' S 1k ) 



Pgmv 



K 



However, if we choose to impose the short sale constraint, i.e., pi > 0, Vi = 1, . . . , K, 
the problem cannot be solved analytically anymore and it requires numerical opti- 
mization techniques. 

The optimal portfolio with (p*, rf) minimizes the risk arising from fluctuations 
in the level of returns of individual assets, i.e. specific risk, generated by spot price 
fluctuations. However, it is still susceptible to the changes in the level of the stock 
market. Therefore, the investor can protect herself by hedging her portfolio. Nor- 
mally a hedge consists of a risky asset and an offset of another related security, 
mostly futures and options. A portfolio can be hedged using futures of a certain 
market index, that portfolio is mirroring (see application using real data in Section 
51). 



As described in Hull (2012), "(...) a futures contract is an agreement between 
two parties to buy or sell an asset at a certain time in the future for a certain price". 
The asset can be either a commodity, or a financial asset, such as currencies of stock 
indices for example. Hedgers constitute a big part of the participants in the futures 
markets, because futures is an effective tool to reduce risk. An agent, who makes a 
deal to sell an asset has a short futures position, therefore it is called short hedge, 



and vice- versa. See Hull (2012) for more on hedging strategies using futures. 

The optimal proportion of the futures contract that counterbalances the spot 
position is called the optimal hedge ratio, D*. Say that rf p is the total return of 
the hedged portfolio at time t, rf is the return of the future contract and rf is the 
return of the portfolio or asset of interest we want to hedge. The total return rf p 
is the difference between the portfolio return and the futures return, scaled by the 
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hedging strategy D: 



MP 



D x ri 



In some cases, commodities for example, the conventional hedging assumes unit 
correlation between the underlying asset and the financial derivative (hedging in- 
strument), where the optimal strategy results in the hedging ratio equal to one. 
Nonetheless, in practice this approach has certain limitations, especially in other 
than commodity markets, see DeCovny and Tacchi (1991). Therefore, the investors 



should employ such a hedging strategy that minimizes the total variance of the 
hedged portfolio: 



D 



GMV 



argminVar [ 



HP~\ 
t J > 



with solution: 



D 



GMV 



Cov [r[,rf] /Var [rf] = p F)P ■ a P /a F . 



However, this approach relies on unconditional variance, and is criticized in a number 
of papers, because the dependence between the futures and the underlying asset is 



likely to be time- varying. Park and Bera (1987) notice that simple OLS approach 
ignores heteroscedasticity. 

The minimum total variance hedging quality criterion is not the only one that 
can be considered in searching for the optimal hedge ratio. Another is a utility-based 



approach, which maximizes hedger's preferences, as seen in Delbaen et al. (2002) 

] 



Becherer (2004), Rossi and Zucca (2002), among others. As before, assume a risk- 

= E [r 



MP 



averse investor, who has a quadratic utility function E [U(rf , 7 
Var [rf p ] . Therefore, solve for the optimal hedge ratio: 



D*u = argmaxE [t/(rf p , 7 )] = argmaxE [r£ p ] - |Var [r£ p ] , 



MPl 



7 /2- 



with solution: 



Dt 



7 ■ Cov(rf,rf) - a 1 
7 ■ Var [rf] 



2.2 Portfolio Allocation and Hedging: Dynamic Approach 

The use of the time-varying covariance matrix to determine portfolio weights and 
hedging ratios leads to better performing portfolios, as shown by 
and reduces the hedged portfolio risk, as seen in 



Choudhry (2004). Giamouridis 



Yilmaz (2011) 



and Vrontos (2007) find that portfolios, constructed under dynamic approach, have 



lower average risk and higher out-of-sample risk-adjusted realized return. 

To solve the portfolio allocation problem in our case, instead of E = Cov[rt] we 
use estimated one-step-ahead conditional covariance matrix for the assets returns 
Cov [r t+ i|X f ] = -f/f+i, which is adjusted continuously on the basis of available infor- 
mation up to time t: X t . Therefore, we are able to obtain optimal portfolio weights 
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for each period: 



Pu,t+i 



\L 



* \T 
VOMV.t+U^t 



1 iV 1 u »' H ^ lK 



7 



■ H t+1 l K 



t+l L K 



l' K H t+1 l K 



(2) 



where fi is a mean vector of the returns, which is assumed to be constant. Same 
goes for subsequent portfolio hedging problems. We can extend our analysis into 
estimating time- varying hedge ratio, that is adjusted at every period: 



7 • Cov [rf ,rf 


L 


]-» F 


L 


7 • Varfrf 


L] 



1 -H 



(2,2) 
t+l 



D 



GMV,t+l 



\L 



Cov [ 



't+l 5 't+l 



It] _ Ht 



(1,2) 
t+l 



Var[rf +1 |X,] H 



(2,2) ' 
t+l 



(3) 
(4) 



where here H t +i is the one-step-ahead covariance matrix between the portfolio we 
want to hedge and the financial instrument, futures in particular. 

In all portfolio and hedging solutions, the H t+ i is one-step-ahead conditional 
covariance matrix of the returns, that is estimated using some multivariate volatility 
model. One possibility is the use of multivariate GARCH models, since they are 
easy to implement and can capture the stylized facts, that are characteristic to 
financial returns. The use of MGARCH models in optimal allocation context was 
first suggested by Cecchetti et al. (1988). Since then, there has been a number 



of papers investigating the differences in estimated hedging ratios and evaluating 
their performance using various approaches, from simple OLS, to bivariate vector 
autorregression (VAR), to GARCH. They show that the use of GARCH-type models 



leads to the overall portfolio risk reduction, see Rossi and Zucca (2002), Kroner and 



Sultan (1993) and Yang and Allen (2004), among others. 



Therefore, given a multivariate time series vector of returns r t : 



r t = n + a t 
H t = Cov[r t |2^ a 



IT 1 / 2 



Hl /2 Cov[e t ](Hl /2 Y, 



where a t is mean-corrected returns, e t is a random vector, such that E[e 4 ] = and 
Cov[e t ] = Ik- There is a wide range of MGARCH models, where most of them differ 
in specifying H t . In the next section we describe the MGARCH model used in this 
paper. 



3 Model, Inference and Prediction 

This section describes the asymmetric multivariate GARCH model used for mod- 
eling volatilities, the implementation of Bayesian non-parametric inference and the 
methodology of obtaining predictive densities of the returns and volatilities. 
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3.1 Asymmetric Generalized DCC Model 



Financial returns exhibit two types of asymmetries: in individual volatilities and 
in conditional correlations. In order to incorporate asymmetric volatility effect, for 



individual time series we choose the GJR-GARCH model by Glosten et al. (1993) 



To model joint volatilities we use Asymmetric Generalized DCC (AGDCC) model, 
proposed by Cappiello et al. (2006) (based on the previous work by |Engle| ( [2002al )). 
We also incorporate the ideas of Hafner and Franses (2009), where the parameters 



in the correlation equation are vectors, not scalars, thus allows for asset-specific 
dynamics. This leads to the following final model: 

l li 

r t = fi + H t ' e t , where r t - \i = a t and e t ~ Fk, (5) 

H t = D t R t D t , (6) 
A 2 = di ag(^) + [diag(ofi) + diag(diag(0 J )/ t , _ 1 )] d t _ x a t _ x + diag(ft) A-i> (?) 

e t = D; l a u 7] t = e t Q I(e t < 0), (8) 

Qt = 5(1 - - A 2 - P/2) + kk' e'^st-x + XX' Q t -x + 55' ^.^-i, (9) 

Rt = (diag(g t ))" 1/2 ^(diag(Q i ))- 1/2 , (10) 

where J-'k is a i^-dimensional distribution, specified later, "diag" stands for either 
taking just the diagonal elements from the matrix, or making a diagonal matrix 
from a vector, S is a sample correlation matrix of e t , denotes Hadamard matrix 
product operator. Parameters k, X and 5 are K x 1 vectors, R 
X = K- 1 Zf =1 A, and 5 = Zf=i 

Q t , we impose Ki, Aj, 5j > and /t 2 + A 2 + 5f/2 < 1, Wi = 1 
volatilities are represented in the equation (J7|: df it = Ui + (ai + 4>ilt-i)af t _ 1 + (3df it _ 1 , 
where / is an indicator function (an < 0) and df it are individual asset volatilities, 
following a GJR-GARCH model with parameters u>i > 0, < ai,<f)i,(3i < 1, such 
that cti + & + <l,Vi = l,...,K. 



K 

i=l K ii 

5i. To ensure positivity and stationarity of 
2 + bill < 1, Vi = 1 K. Individual 



The AGDCC by Cappiello et al. (2006) is just a special case where K\ 
kk, Ai = . . 



Xk and 5% 



Qt = 5(1 



K 



X - 5/2) + kx e' t _ x e t -\ + A x Q t _ x + 5 x r}' t _ 1 r] t -x. (11) 



As for the distribution of e t ~ Fk, we model it as an infinite scale mixture of 
Normal distributions, where the density of e% is as follows: 

h{e t \G)= /A^(e i |0,A i )rfG(A- 1 ), 



where A/x(q|0, A^ 1 ) denotes a f^-dimensional conditional density function of multi- 
variate Normal distribution with mean zero and scale matrix A t , and G is the scale 
mixing distribution, which is unknown and modeled by a Dirichlet Process. 

DP is a multi-parameter generalization of the Beta distribution and defines a 
distribution over distributions. DP leads to discrete probability measures, which is a 
disadvantage in practice. This problem can be overcome by using Dirichlet Process 
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Mixture model, as seen in Antoniak (1974): 



e t |A t -A^^A- 1 ), 



A.AG ~ G, 



(12) 



G\c, G ~DP(c,G ) 
c ~ vr(c), 



where c > is a scale parameter with prior density ir, and Go is a prior probability 
measure. Observe that G is a random distribution drawn from the DP and because 
it is discrete, multiple A,'s can take the same value simultaneously, making it a 
mixture model. Hence, using the stick-breaking representation, the hierarchical 



model in (12) can be seen as an infinite mixture of distributions: 



f(e t \A,p) = ^p l AV(e t |0,A: 



(13) 



where the weights are obtained as follows: p\ = vi, pi = (1— v%) . . . (1— Vi-i)Vi, where 

Vi is Beta distributed: v , ~ B (1, c) for i — 2, 3, We assume a conjugate model, 

where Go is a Gamma distribution with parameters (a/2,6/2) and also assume a 
Gamma hyper-prior on the concentration parameter c ~ Q{do, bo). Finally, we have 
assumed uniform prior distributions for the parameters of the GJR-AGDCC model. 



3.2 Bayesian Inference 



A number of papers in the field of GARCH-type models have explored different 
Bayesian procedures for inference and prediction and different approaches to mod- 
eling the fat-tailed errors and/or asymmetric volatility. The recent development of 
modern Bayesian computational methods, based on Monte Carlo approximations, 
such as importance sampling, and MCMC methods, such as the Metropolis-Hastings 



algorithms, have facilitated the usage of Bayesian techniques, see Robert and Casella 



(2004). 



The following section describes the Bayesian non-parametric procedure used for 
the GJ R-AGDCC model in (|5])-( 10). The algorithm is based on works by Walker 



et al. (2011). 



(2007), Papaspiliopoulos and Roberts (2008), Papaspiliopoulos (2008) and Ausin 



Regarding the inference algorithms, there are two main types of approaches. On 
the one hand, the marginal methods, which rely on the Polya urn representation. 
All these algorithms are based on integrating out the infinite dimensional part of the 
model. One of the most recent papers, based on this method in MGARCH setting 
is by Jensen and Maheu (2012). Recently, another class of algorithms, called condi- 



tional methods, have been proposed. These approaches, based on the stick-breaking 
scheme, leave the infinite part in the model and sample a finite number of vari- 



ables. These include the procedure by Walker (2007), who introduces slice sampling 



schemes to deal with the infiniteness in DPM, and the retrospective MCMC method 



of Papaspiliopoulos and Roberts (2008), that is later combined by Papaspiliopoulos 



(2008) with a slice sampling method by Walker (2007) to obtain a new composite 
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algorithm, which is better, faster and easier to implement. Generally, the stick- 
breaking procedure, compared to the Polya urn, produces better mixing and simpler 
algorithms. 

As seen in Walker (2007), in order not to sample an infinite number of values 
at each MCMC step, we introduce a latent variable Ut, such that the joint density 
of (e, u) given (p, A) is given by 

OO 

f(e t ,u t \p,A) = < P i)N K {et\^K 1 )- ( 14 ) 

i=i 

Let A p (ut) = {i : pi > Ut} be a set of size N Ut , which is finite for all u% > 0. Then 
the joint density of (et,Ut) in (14) can be equivalently written as f(et, Ut\p, A) = 
S«eA (u t ) -A/}^(e t |0, Aj -1 )- Integrating over u t gives us the previous density of infinite 
mixture of distributions (13). Finally, given u t , the number of mixture components 
is finite. In oder to simplify the likelihood, we also need to introduce further indi- 
cator latent variable z t , which indicates the mixture component that e t comes from: 
f(e t ,z t — j, u t \p, A) = A/V(et|0, A _1 )l(j G A p (u t )). Define parameter sets Q = (p,A) 
and $ = (p, oj, a, (3, <fr, n, A, 5), where = (Q, $) is the set of all model parameters. 
Then, the log likelihood of 0, given the latent variables u t and z t looks as follows: 

1 T 

l(e\ut,zt) = --£>log(27r) + log \H;\ +a t H* t - 1 a t ), (15) 
t=i 

where H£ is the new conditional covariance matrix, adjusted by the variance of the 
errors: 

Cov folZU, ^) = H t 1/2 A;X /2 = H l 

Next, we describe the DPM model algorithm step by step. 

Firstly, given z, the conditional posterior distribution of concentration param- 
eter c is independent of the rest of the parameters, as in Escobar and West (1995). 
So, we first sample an auxiliary variable £ ~ £>(c + 1, T) and then c from a Gamma 
mixture: 

7r^(a + z*,b - log(O) + (1 - ^)£(a + z*-l,b - log(O), 

where z* = max(zi, . . . ,zt) and = (a + z* — l)/(a + z* — 1 + T(b — log(£))). 

In the second step sample the weights of the components Vj for j = 1, . . . , z* , 
where the prior for v ~ B(l, c) and, given the data and z: 



Vj\Z ~ 



B(nj + l,T-^2m + c), 

i=i 

j 

where rij is the number of observations in the j th component and ^2 n\ gives the 

l=i 

cumulative sum of the groups. Also, p\ = V\, pj = (1 — Vi) ... (1 — Vj-i)Vj, for 
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At the third step update ut ~ U(0,p Zt ), for t = 1, . . . ,T. 

In the fourth step sample all the values of pj that are larger than u t . As 



Walker (2007) showed we need to find the smallest j* such that Yjj=iPj > u * an d 
then update Vj and pj for j = z* + 1, . . . , j*, where u* = min(wi, . . . ,Ut)- 

In multivariate setting the conjugate prior for the inverse of covariance matrix 
is the Wishart distribution, which can be seen as a matrix generalization of the 



chi-square distribution (Eaton, 2007). Therefore, the prior for the scale matrix A is 



chosen to be Wishart (V,df,K), where V is positive semi-definite symmetric scale 
matrix of dimensions K x K, df is degrees of freedom and K is the dimension 
parameter. The real degrees of freedom can be obtained as df + K — 1. The mean 
vector of the variables distributed as Wishart is E [A] = (df + K — 1) x V. Sampling 
variability is large when the degrees of freedom is small. Thus, df has to be larger 



than K—l, and V such that the expectation of the At is identity: V 



i 



seen m 



(df+K-l) 



Ik, as 



Jensen and Maheu (2012). See more on the properties of Wishart distribution 



in Eaton (2007). Next, update the A, whose posterior distribution is independent 



of (p, u t ) 



A,|- ~ 



W (df + n j ,(V~ 1 + Cov[e Zt=j \) 



In the sixth step update to which component the observations belong to by 
using the following (as seen in |Walker| ( |2007l )): 



Probability^ = cx 1 (j G A p (u t )) A/k(e t |0, A, 1 ), 



where A p (ut) = {j : Pj > Ut}, which is not empty. 

The rest of the steps of the algorithm concern updating the parameters of the 
GJR-AGDCC model. We use the Random Walk Metropolis Hasting (RWMH), 
where for each parameter 9 G a candidate value 9 is generated from a A'-variate 
Normal distribution with mean equal to the previous value of the parameter and 
variance calibrated to achieve a desired acceptance probability. This procedure is 
repeated at each MCMC iteration. The probability of accepting a proposed value 

9, given the current value 9, is a(9,9) = min jl, n^=i K r t\9) / K r t 1^)}; where 
the likelihood used is as in (15), see e.g. Robert and Casella ( 2004[ ). In this paper 
the acceptance probabilities are adjusted to be between 20% and 50%. For more 



detailed explanation of the algorithm in univariate setting see Ausin et al. (2011). 



3.3 Prediction 

In this section, we are mainly interested in estimating the one-step-ahead predictive 
density of the returns: 



f(r t+1 \n, ...,r t ) 



f(r t+1 \^)f(e\r 1 ,...,r t )dQ. 



(16) 



Although this integral is not analytically tractable, we can approximate it using the 
MCMC output. For this, we make use of the the procedure described in Walker 
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(2007), where it is explained how to obtain a sample from the predictive density of 
the errors f(et+i\e\, . . . , e<). Thus, having this sample and a posterior sample of the 
remaining model parameters $, it is easy to obtain a sample from the predictive ( 16 ) 
using equations (|5])-( 10 ). The following is the detailed explanation of the procedure. 

At each MCMC iteration there are weights pj and corresponding precision ma- 
trices Aj. We sample a random variable r ~ U(0, 1) and take such Aj for which 
Pj-x < r < pj. If we need more weights (i.e. pj < 1), we can sample additional pj 
as before^J and Aj from the prior. Once we have Aj selected, Walker (2007) proposes 
to draw one observation from each Af (0, Aj 1 ) in order to have a sample {et+i,™,}^ 1 



-il e i) 



1 m=l 
■6*)- 



of the size of the MCMC chain from the predictive distribution of f(e u 
Alternatively, instead of sampling just one observation, we suggest to assemble the 
entire collection of precision matrices, {A Zt+i m } . This approach has several ad- 
vantages. Firstly, it allows to incorporate prior information about the variance of 
the errors by having the small probability of being sampled from the prior. Secondly, 
it allows to increase the sample size for the predictive density of /(e t+ i|ei, . . . , e t ), 
because instead of sampling just one observation as in |Walker| ( |2007l ), we can sample 

as many as we choose at each MCMC iteration: q ~ Af(0, A -1 ), % = 1, . . .. And 
finally, it provides a sample of one-step-ahead volatilities {H£ +1 m } M ., that we will 
use in the portfolio allocation and hedging problems in the following section. 



3.4 Bayesian Portfolio Decisions 

As commented in the introduction, optimal allocation is greatly affected by the pa- 



rameter uncertainty, which has been recognized in a number of papers, see Jorion 
(1986) and Greyserman et al. (2006), among others. They conclude that in fre- 



quentist setting the estimated parameter values are considered to be the true ones, 
therefore, the optimal portfolio weights tend to inherit this estimation error. Instead 
of solving the optimization problem on the basis of the choice of unique parame- 
ter values, the investor can choose the Bayesian approach, because it accounts for 



parameter uncertainty, as seen in Kang (2011) and Jacquier and Poison (2012), for 
example. 

Portfolio decision problems in Bayesian setting are usually solved by choosing 
such portfolio weights, that maximize the expected utility of the portfolio with 
respect to the predictive density of the one-step-ahead returns: 



4 = argmax 

p 



U(p'r t+1 )f(r t+ i\r 1 , . . . ,r t )dr u 



Therefore, the investor would obtain point optimal portfolio weights, where the 
parameter uncertainty has been accounted for. In our case, we do not have the 
analytically tractable posterior distribution of the returns, just a sample of size M. 



Pi 



v i, Pj = (l-v 1 )...{l 



8(1, 
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We can approximate the solution: 



M 



Pt+i ~ argmax ^- VJ7(pV t+ i im ), 

m=l 

where {rt +l m }^ =1 is a predictive sample of one-step-ahead returns, obtained as 
explained in the previous section. However, this approach provides only with a 
point estimation of the optimal portfolio weights. Since the analytically tractable 
posterior distribution is not available, it is not straightforward to obtain measures of 
uncertainty as credible intervals in order to asses the quality of this estimation, see 



Brandt (2009), for example. Moreover, this approach does not provide the measure 
of uncertainty for subsequent portfolio characteristics, such as the portfolio expected 
return, variance or expected utility. 

Alternatively, we propose to obtain samples from the entire posterior dis- 
tribution of optimal portfolio weights /(p* + i|ri, . . . ,r t ). This approach relies 
on solving the allocation problem at every MCMC iteration and approximate for 
example the posterior mean of the optimal portfolio weights by: 



f 1 M 

E\p* +1 \ ri , ...,r t }= / p* +1 f{e\ ri , . . .,r t )de « — 2^P t * +1 , m . 

m=l 



where {Pt+im}m=i ls a posterior sample of optimal portfolio weights obtained for 
each value of the model parameters in the MCMC sample. In other words, since 
we have assembled M one-step-ahead volatility matrices and mean vectors, we can 
solve the portfolio allocation problem M times. Similarly, we can approximate 
the posterior median of p* t+l and credible intervals by using the quantiles of the 
sample of optimal portfolio weights. In this manner, we are able to obtain a sam- 
ple from the posterior distribution of portfolio expected returns {E[r^_i]m}m=i = 
{Pt+i, m /i m }*=i, variance {Var[r^ 1 ] m }* f =1 = {K+i, m ^+i,mK+i, m }m=i and expected 
utility {E[U t+1 ] m }% =1 . 

In hedging exercise, same as in the portfolio allocation problem above, it is 
possible to obtain a sample from the posterior distribution of the optimal hedge 
ratios D* +1 . As for the point estimate, the investor can choose either the posterior 
sample mean or median, which can be obtained from the collection of optimal hedge 
ratios {D* +1 }% =v 



4 Simulation Study 

The goals of this simulation study are to show the flexibility and adaptability of the 
DPM model and to explain some bimodal posterior densities that we later observe 



in real data applications. Here we use the basic GJR-AGDCC model as in (11) 



by Cappiello et al. (2006). We have generated three bivariate time series of 3000 
observations with the following errors: Gaussian Af(0,I 2 ), Student-t T{I 2) v = 8) 
and a mixture of two bivariate Normals 0.9jV(0, a% = 0.8, oyi = 0.0849, o 2 = 0.9) + 
0.1jV(0, <j\ = 2.8, <7i2 = —0.7637, a 2 = 1.9). In the mixture data we have chosen 
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a bigger variance of the second component for the first series, to make it more 
volatile than the second. Then, estimate all three data sets using a DPM model. 
The point estimates are not reported in the paper because of the limited space. All 
parameters were estimated well, with true parameters always inside the 95% credible 
intervals. The length of MCMC chain is 10,000 burn-in plus 20,000 iterations. The 
contour plots in Figure [T] compare true predictive densities of returns f(r t +i\It) 
with the estimated ones, which were obtained by sampling 5 observations at each 
MCMC step, resulting in a sample size of 100,000, as explained in the previous 
section. As we can see, the estimation results are very precise compared to the true 
contour of one-step-ahead returns. The contours can be seen as a summary of the 
estimation results for all 13 parameters of the model $ = (/i, u, a, f3, <p, k, A, 5) and 
the distribution for the error term. Therefore, the infinite mixture model is a flexible 
tool that is able to adjust to whatever distribution the data comes from. 

Figure [T] goes here 

Next, the top part of the Figure [2] presents the marginal kernel smoothing 
densities for one-step-ahead errors, f(et+i\ri, . . . , r r ). The densities are symmetrical, 
thus we present only the left and right tails for each marginal series. The Student-t 
data predicts fatter right and left tails, consequently, allowing for more extreme 
observations, which in turn increase the volatility. The mixture model predicts fat 
tails just for the first marginal series, but not for the second. This is due to the 
nature of the simulated data: the second mixture component is more volatile for 
the first series than for the second (a\ = 2.8, <Ji = 1.9). The effect of how fat 
tails increase volatility can be seen in the middle row of Figure [2j where the kernel 
smoothing densities of the elements of the error covariance matrices {^-t+im\m=i 
are presented. The Gaussian error data does not allow for extreme observations, 
therefore does not allow for high volatilities. However, in a Student-t and mixture 
data the predictive density of the variance of the errors has a very fat right tail, 
which means high volatilities. Because of the more extreme returns, the density of 
the variance of the errors is not symmetrical. And finally, the bottom row of Figure [2] 
presents the predictive densities of the volatilities of the returns {H^ +1 m }^f =1 . Here 
the effect of the fat tails is even more pronounced, since we also take into account 
the asymmetric volatility effect. The middle graph, which in the middle row was 
symmetrical, here is also bimodal, because it models the asymmetric correlation 
effect. This simulation study helps to understand and explain the bimodal posterior 
distributions that appear in the real data application in the next Section. 

Figure [2] goes here 

5 Data and Results 



In the illustration using real data, same as in the simulation study, for the sake of 
simplicity we use the basic GJR-AGDCC model, where k, A and 5 are scalars as in 
(fTTl) bylCappiello et al\ (120061). 
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In this section we illustrate the financial applications described in Section |2j 
At first we solve portfolio allocation problem for bivariate time series for utility- 
based and GMV approaches. Then, once we have the portfolio, we find the optimal 
hedge ratios under minimum variance and maximum utility criteria using futures of 
a certain index. All data was obtained from IHS Global Insight database and Yahoo 
Finance. 

5.1 Portfolio Allocation 

For the first illustration we use the daily price data of Apple Inc. company (P^) and 
NASDAQ Industrial index (P t N ) from January 1, 2000 till May 7, 2012. Then, daily 
prices are transformed into daily logarithmic returns (in %), resulting in 3098 ob- 
servations. Table [T] provides the basic descriptive statistics, and Figure [3] illustrates 
the dynamics of returns. 

Table 1. Descriptive Statistics of the Apple Inc. and NASDAQ Ind. Return Series 

100 x In JJpj 100 x In ^§Tj 



Mean 0.0973 0.0020 

Median 0.1007 0.0766 

Variance 9.7482 3.1537 

Skewness -4.2492 -0.1487 

Kurtosis 102.0411 7.1513 
Correlation 0.5376 



Figure [3] goes here 

As expected, the Apple Inc. has higher overall variance because of the higher 
mean return. Both returns do not exhibit any evidence of auto-regressive behavior. 
Apple Inc. returns contain one atypical data point, corresponding to September 29, 
2000. The very low return is due to an announcement the day before about lower 
than expected sales. The estimation results of the AGDCC model are reported 
in the Table [2j The ML and RWMH estimation approaches, assuming Gaussian 
innovations, as expected, provide very similar estimates. DPM model exhibits dif- 
ferences in point estimations for the parameter, which measures the asymmetric 
volatility effect for marginal return series. Also, u parameter for both Gaussian 
models is larger than in DPM, probably because of the atypical data in the return 
series. Gaussian error models have to inflate the constant volatility parameter in 
order to accommodate this atypical data point, whereas in the DPM model is has 
been accounted for in one of the mixture components of the distribution of the 
errors. 
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Table 2. Estimation Results for Apple Inc. (1) and NASDAQ Ind. (2) Returns, 
20,000 iterations plus 10,000 burn-in 





ML-Gaussian Errors 


Bayesian Gaussian Errors 


Bayesian DPM 




Parameter 


Post, mean 


Post, mean 




(st.dev.) 


(Post. st. dev.) 


(Post. st. dev.) 




0.0973, 0.0020 


0.1628,0.0231 


0.1404, 0.0417 




(0.0296), (0.0149) 


(0.0408), (0.0216) 


(0.0362), (0.0215) 




0.1751, 0.0232 


0.2676, 0.0277 


0.1324, 0.0206 




(0.0142), (0.0024) 


(0.0599), (0.0053) 


(0.0450), (0.0055) 


ax, ot-2 


0.0672, 0.0059 


0.0923, 0.0146 


0.0692, 0.0125 




(0.0050), (0.0047) 


(0.0135), (0.0073) 


(0.0149), (0.0066) 




0.8725, 0.9250 


0.8396, 0.9226 


0.8893, 0.9250 




(0.0042), (0.0053) 


(0.0166), (0.0081) 


(0.0173), (0.0078) 


01,02 


0.1095, 0.1189 


0.1090, 0.1016 


0.0506, 0.0804 




(0.0094), (0.0075) 


(0.0252), (0.0124) 


(0.0234), (0.0168) 


K 


0.0231 


0.0113 


0.0140 




(0.0027) 


(0.0090) 


(0.0078) 


A 


0.9628 


0.9794 


0.9516 




(0.0047) 


(0.0195) 


(0.0314) 




0.0061 


0.0051 


0.0182 




(0.0049) 


(0.0042) 


(0.0117) 



After the estimation was carried out using the Gaussian errors and the infinite 
mixture of Gaussian distributions, we are able to approximate the predictive one- 
step-ahead return distributions f [r t+ i\X t ). Figure [4] shows the tails of the kernel 
smoothing marginal densities of the one-step-ahead returns. We can observe the 
differences in tails arising from different specification of the errors. The DPM model 
permits for a more flexible distribution, therefore, for more extreme returns, i.e. 
fatter tails. The kernel smoothing densities were obtained using the same procedure 
as in the simulation study: by sampling 5 observations at each MCMC step, resulting 
in a sample size of 100,000. 

Figure [4] goes here 

Table [3] presents the estimated mean, median and 95% credible intervals of one- 
step-ahead volatility matrices in Bayesian context, obtained from the collection of 
M matrices {H* +1 m }^f =1 . Both Gaussian models provide very similar volatility esti- 
mates. In the constant unconditional approach, the individual volatilities and covari- 
ances are much greater. Unconditional covariance matrix gives estimated volatility 
of 9.75 for Apple Inc. and 3.15 for NASDAQ Ind. with the correlation of 0.5375, 
whilst the dynamic Gaussian and DPM models for t + 1 estimate 7.42 and 7.50 for 
Apple Inc. and twice as small for NASDAQ Ind. 1.59 and 1.56 correspondingly, with 
correlations of 0.5329 and 0.4972 respectively. In Gaussian model the means of the 
marginal volatilities are very similar to the medians, meaning, that the distribu- 
tions are symmetric. In DPM model the means are greater than the medians, which 
means that the posterior distributions are skewed to the right, i.e. they have longer 
right tails. This happens because DPM model permits some returns to come from a 
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very volatile component, as seen in Figure [4] and explained in the simulation study 
in Section m Fi gure p\ illustrates the right tail of the posterior distribution of one- 
step-ahead volatilities of Apple Inc. returns, where, P (h*^ 1 ^ > 13.45 \Xt\ = 0.050, 

P \ H*+-y^ > 60.51 \%t\ = 0.010. This reads, that there is 1% chance to observe a 

data point, that is more volatile than 60.51 and 5% chance to observe volatility 
greater than 13.45. 

Figure [5] goes here 



Table 3. Estimated Means, Medians and 95% Credible Intervals of One-Step- Ahead 
Volatilities of the Apple Inc. and NASDAQ Ind. Return Series 



Matrix element 




(1,1) 


(1,2), (2,1) 


(2,2) 


Constant 


E 


9.7482 


2.9805 


3.1537 


ML Gaussian 


Ht+i 


7.4631 


1.6646 


1.5729 


Bayesian Gaussian 


mean 


7.4185 


1.8330 


1.5947 




median 


7.4272 


1.8225 


1.5930 




95% 


6.9514, 7.8445 


1.6572, 2.0515 


1.4759, 1.7237 


Bayesian DPM 


mean 


7.5016 


1.7010 


1.5605 




median 


4.7396 


1.5844 


1.4810 




95% 


3.2745, 44.9882 


0.8202, 3.4372 


0.9724, 3.1264 



Utility and GMV portfolios. Here we solve for the utility-based and GMV 
portfolios without the short-sale constraint, as in Q-fl). The risk 

aversion coef- 
ficient 7 = 0.03, which means that the penalty for increased variance is 7/2. The 
choice of the risk-aversion coefficient is arbitrary. Table [4] shows the estimated means 
and credible intervals of optimal portfolio weights p$ +1 , portfolio return r[ +l , port- 
folio variance ofi^ and investor's utility from the portfolio U[ +1 . In the utility-based 
case, the optimal portfolio weights for constant model are very different for the rest 
of the models: it suggests investing 48% of the wealth in Apple Inc. shares, and 
the rest in the index. Gaussian and DPM models estimate more similar portfolio 
weights, 0.83 and 0.94, where both mean estimates enter in the credible intervals 
of the others. Gaussian 95% credible interval is around two times smaller than the 
DPM, which is expected, since DPM model permits for fatter tails in returns and 
volatilities. The rest of the columns present the estimated expected returns, vari- 
ances and utilities for the utility-based portfolio. Gaussian and DPM models result 
in very similar point estimates, whilst the DPM model permits for fatter tails in all 
cases. The 95% credible intervals for portfolio returns, variances and utility is from 
25 to 30% wider in the DPM setting. This illustrates, that if the investor chooses 
to impose a restrictive distribution of the errors, she would not be able to measure 
the uncertainty of her financial decisions correctly. In Gaussian setting, she would 
be overconfident. Similar conclusions can be drawn from the GMV-based portfolio 
estimation results. The constant covariance model proposes investing 2.5% of the 
funds in Apple Inc. shares, however, the total portfolio risk is as twice as big as 
in Gaussian and DPM models. Both models propose short-selling the risky asset 
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and investing all the funds plus the income from short-selling into a less risky asset 
- the index. In GMV-based portfolio the differences in uncertainty are even more 
obvious: the credible interval for p*GMVt+\ f° r DPM model is more than three times 
wider. The general differences in tails can be observed in the Figures [6] and [7j 



Table 4. Estimated Utility-Based and GMV Portfolio Weights of the Apple Inc. 
and NASDAQ Ind. Return Series 







Pu,t+i 






P*GMV,t+l 






Constant 


Bayesian Gauss 


Bayesian DPM 


Constant 


Bayesian Gauss 


Bayesian DPM 


* 

Pt+1 


0.4825 


0.8253 


0.9359 


0.0249 


-0.0453 


-0.0312 


(95%) 




(0.4070, 1.2277) 


(0.0682, 1.9681) 




(-0.0818, -0.0132) 


(-0.1833, 0.1215) 


p 
r t+i 


0.0479 


0.1456 


0.1425 


0.0044 


0.0168 


0.0386 


(95%) 




( 0.0406, 0.2807) 


(0.0324, 0.3332) 




(-0.0258, 0.0604) 


(-0.0055, 0.0828) 


(T 2P 


4.6022 


5.8757 


4.9925 


3.1493 


1.5822 


1.5301 


(95%) 




(2.6560, 10.2364) 


(1.7159, 11.1193) 




( 1.4600, 1.7131) 


(0.9398, 3.0042) 


u t+i 


-0.0211 


0.0575 


0.0676 








(95%) 




(-0.0074, 0.1302) 


(-0.0095, 0.1704) 









Figures [6] and [7] present the kernel smoothing densities of portfolio returns, 
variances and utilities for Gaussian and DPM models. In all cases DPM model 
exhibits fatter right tail and wider credible intervals, which arises as a consequence 
of the fat-tailed return distribution. 

Figure [6] goes here 

Figure [7] goes here 

To sum up, these portfolio allocation exercises helped to illustrate the direct 
consequences of return distribution to the uncertainty of financial decisions. DPM 
model permits the investor to perform inference and prediction about the returns 
and their volatilities without imposing arbitrary restrictions on the data generating 
process. 



5.2 Hedging 



Once the optimal portfolio has been selected, the investor might want to eliminate 
risk associated from the movements in the market. Say, that the S&P 500 Index 
represents our well-diversified portfolio, therefore, S&P 500 Index futures can be 



used to hedge it. As Hull (2012) notices, hedging using index futures eliminates 



the market risk and the performance of portfolio depends only on the performance 
relative to the market. 

For this illustration we use the daily price data of S&P 500 index (P/) and E- 
mim S&P 50(Q (from now on, just E-mini) (P/) from January 1, 2000 till May 11, 

2 E-mini S&P 500 is a stock market index futures contract, traded on the Chicago Mercantile 
Exchange's Globex electronic trading platform. The value of the future contract is 50 US dollars 
times the E-mini S&P 500 futures price. It is a smaller version (l/5 tft ) of regular S&P 500 futures 
contract, made to be available to smaller investors. 
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2012. The E-mini prices are quoted for the second expiration future nearby. Daily 
prices are transformed into daily logarithmic returns (in %), resulting in 3109 ob- 
servations. Table [5] provides the basic descriptive statistics, and Figure [8] illustrates 
the dynamics of returns. 

Table 5. Descriptive Statistics of the S&P 500 Index and E-mini Futures Return 
Series 



100 x In j£- 100 X In -9 



Mean -0.0063 -0.0023 

Median 0.0592 0.0551 

Variance 1.9417 1.8805 

Skewness 0.0047 -0.1576 

Kurtosis 14.5295 10.1901 

Correlation 0.9792 



Figure [8] goes here 

The estimation was carried out using Gaussian distribution and the infinite 
mixture of Gaussian distributions for the errors. Estimation results are presented 
in Table [6} Both Bayesian approaches provide similar estimates, whereas ML has 
problems in estimating a and 5. Next, Figure [9] presents the predictive kernel 
smoothing densities for marginal returns, obtained by sampling 5 returns at each 
MCMC step, resulting in a sample of 100,000. Observe, that the DPM model allows 
for fatter tails - more extreme returns. 

Table 6. Estimation Results for S&P 500 (1) and E-mini (2) Returns, 20,000 
iterations plus 10,000 burn-in 





ML-Gaussian Errors 


Bayesian Gaussian Errors 


Bayesian DPM 




Parameter 


Post, mean 


Post, mean 




(st.dev.) 


(Post. st. dev.) 


(Post. st. dev.) 




-0.0063, -0.0023 


-0.0159, -0.0142 


0.0222, 0.0219 




(0.0104), (0.0105) 


(0.0150), (0.0153) 


(0.0152), (0.0157) 




0.0188, 0.0166 


0.0199, 0.0206 


0.0161, 0.0174 




(0.0009), (0.0010) 


(0.0023), (0.0026) 


(0.0032), (0.0036) 




3.9 x 10~ 9 , 5.9 x 10~ 9 


0.0126, 0.0208 


0.0165, 0.0224 




(4.6 x 10- 6 ), (1.7 x 10- 6 ) 


(0.0060), (0.0073) 


(0.0069), (0.0082) 


Pi, ft 


0.9059, 0.9142 


0.9099, 0.9046 


0.9122, 0.9084 




(0.0026), (0.0027) 


(0.0068), (0.0082) 


(0.0087), (0.0098) 




0.1583, 0.1458 


0.1162, 0.1123 


0.1076, 0.1075 




(0.0055), (0.0055) 


(0.0094), (0.0096) 


(0.0180), (0.0177) 


K 


0.1336 


0.1402 


0.1153 




(0.0061) 


(0.0153) 


(0.0184) 


X 


0.5802 


0.4114 


0.5276 




(0.0328) 


(0.1067) 


(0.0866) 


s 


4.6 x 10~ 8 


0.0083 


0.0106 




(1.8 x 10~ 5 ) 


(0.0077) 


(0.0100) 
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Figure [9] goes here 

Then, as in the portfolio optimization problem, we are able to obtain one- 
step-ahead volatilities, where the estimated variances and covariances are presented 
in the Table The mean, median and 95% credible intervals of one-step-ahead 
volatility matrices in Bayesian context are obtained from the collection of M ma- 
trices {Ht+im}m=i- Here, the constant unconditional approach estimates variances 
that are twice as big than in dynamic models, for both time series. The mean cor- 
relation is similar in all three cases: 0.9792 for unconditional, 0.9723 for Gaussian 
and 0.9757 for DPM model. The means and medians in Gaussian model are almost 
the same, therefore, the posterior one-step-ahead distribution of volatilities is sym- 
metric. However, in DPM model the means are greater than the medians, meaning 
that the posterior distribution is positively skewed. In order to investigate in more 
detail the posterior of one-step-ahead volatility matrix in DPM model we draw the 



kernel smoothing densities, which are presented in Figure 10 Obviously, they are 
bimodal. As seen in simulation study in Section |4j this can be explained by the 
fact that there are some returns coming from a more than usual volatile component 
(-s). In addition, the asymmetric volatility pronounces this effect even more for 
both marginal variances and asymmetric correlation effect causes the right graph of 



Figure 10 to have bimodal-type form. 



Figure 10 goes here 



Table 7. Estimated Means, Medians and 95% Credible Intervals of One-Step- Ahead 
Volatilities of the S&P 500 and E-mini Return Series 



Matrix element 




(1,1) 


(1,2), (2,1) 


(2,2) 


Constant 


£ 


1.9417 


1.8712 


1.8805 


ML Gaussian 


Ht+i 


0.8851 


0.8128 


0.7885 


Bayesian Gaussian 


mean 


0.8197 


0.7532 


0.7315 




median 


0.8196 


0.7531 


0.7311 




95% 


0.7794, 0.8607 


0.7170, 0.7908 


0.6949, 0.7695 


Bayesian DPM 


mean 


0.8501 


0.7853 


0.7620 




median 


0.7034 


0.6593 


0.6534 




95% 


0.3540, 2.0363 


0.3378, 1.7537 


0.3432, 1.6121 



Utility and GMV Hedging Ratios. Finally, using the estimated one-step- 
ahead volatilities and means, we solve the hedging problem in two cases: using the 
minimum variance criteria and the maximization of expected utility, as in ([3])-(|4]). 
Here we use 7 = 0.3, which means that the penalization for increased variance is 
greater than in portfolio allocation problem, because our principle objective now is 
the reduction of variance rather than gain from the return. Table [8] presents the 
means and 95% credible intervals of optimal hedge ratios D^ +1 , the total hedged 
portfolio return rE^, the hedged portfolio variance and utility U^. The 

mean optimal utility-based and GMV hedge ratios for Gaussian model enters the 
DPM credible intervals, but not vice versa. The medians for utility-based portfo- 
lio are 1.09677 and 0.9194 for Gaussian and DPM models respectively, indicating, 
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that the posterior distribution of optimal hedge ratios for DPM model is positively 
skewed, meanwhile the Gaussian model provides symmetrical posterior densities. 
The same is observed for the GMV portfolio, where the medians are 1.0297 and 
1.0053 for Gaussian and DPM models respectively. Even though the hedging ratios 
are different for both approaches, they result into statistically equal hedged portfolio 
variances and utilities. The constant covariance matrix hedging approach, compared 
to both time-varying Bayesian models, provides worse hedged portfolio with high 
variance and small utility. In all the cases, credible intervals for DPM model are 



from 20% up to six times wider. This also can be observed from Figures 11 and 12 



Table 8. Estimated Utility-Based and GMV Optimal Hedging Ratios of the S&P 
500 and E-mini Return Series 



u u,t+i 



D, 



GMV.t+l 



Constant Bayesian Gauss Bayesian DPM Constant Bayesian Gauss Bayesian DPM 



n* 

(95%) 

T 
r t+l 

(95%) 

2T 
a t+l 

(95%) 

uf +1 

(95%) 



0.9992 
-0.0039 
0.0799 
-0.0159 



1.0947 
(0.9539, 1.2265) 

0.0007 
(-0.0071, 0.0110) 

0.0506 
(0.0423, 0.0723) 

-0.0069 
(-0.0138, 0.0008) 



0.8974 
(0.6089, 1.1006) 

0.0039 
(-0.0040, 0.0180) 

0.0522 
(0.0239, 0.1374) 

-0.0039 
(-0.0193, 0.0064) 



0.9950 
-0.0039 
0.0798 



1.0298 
(1.0178, 1.0422) 

-0.0013 
(-0.0077, 0.0050) 

0.0440 
( 0.0413, 0.0469) 



1.0164 
(0.9479, 1.1] 
-3xl0" 5 
(-0.0065, 0.0062) 

0.0390 
(0.0205, 0.1283) 



Next, Figures 11 and [12] present kernel smoothing densities of optimal hedge 
ratios, hedged portfolio expected return, variance and expected utility for both 
hedging approaches. As commented before, the DPM credible intervals are wider 



due to the fatter tails in the return distribution. In Figure 12 the second and third 
graphs demonstrate strong bimodality of the posterior densities, which is a direct 
consequence of the bimodality in predictive volatility densities. 



Figure 11 goes here 



Figure [12] goes here 

Here we present a short numerical illustration for the portfolio hedging problem. 



As seen in Hull (2012), the optimal number of future contracts is given by N* 



D*(Vp/Vp), where D* is the optimal hedge ratio, Vp is the value of the portfolio 
and Vp is the value of one futures contract. In time- varying setting, the value of the 
portfolio or futures contract depends on the prices at t + 1, that can be obtained 
as P t+ \ = Pt ■ exp{r t+1 /100}, where P t is the current price. Therefore, the optimal 
number of contracts at each MCMC iteration is: 



N* — D* 



shares Pf exp{rf +1 /100} 
^0~'pfexp{rf +1 /100}' 



Say, an agent owns 1,000 shares, the price of one share of the portfolio is 1, 412.00 
and the E-mini futures price for the second expiration contract nearby at the mo- 
ment is 1, 400.00. The optimal number of contracts for constant covariance model is 
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20.0715, Gaussian - 20.7727, DPM - 20.5022. In Bayesian models the optimal num- 
ber of contracts is a sample mean of MCMC chain, whose posterior kernel smoothing 



densities can be seen in Figure 13 So, if the investor employs a time-invariant ap- 
proach, she would short A^* +1 = 20 E-mini futures contracts, however, in order to 
minimize the portfolio variance at t + 1, the investor should short iV t * +1 = 21 in 
Gaussian and N£ +l = 21 in DPM. The 95% credible interval for Gaussian model is 



(20.5298, 21.0221) and for DPM - (19.1196, 22.5682), which is 7 times wider. Since 
Gaussian errors do not permit extreme returns and high volatilities, the Gaussian 
credible interval for N* +1 is smaller than in reality, making the investor overconfident 
about her decisions and ignore risk arising from the extreme returns. 



Figure 13 goes here 



In portfolio allocation or in hedging context, adjusting portfolio weights at each 
period might lead to high transaction costs, thus the investor will adjust her portfolio 
only if the expected utility after the adjustment minus the transaction costs is greater 
than the expected utility without the adjustment. Both illustrations have shown 
the differences in error specifications in using real data. We have illustrated how 
quantification of uncertainty reflects distributional assumptions of the errors. 



6 Conclusion 

In this paper we have considered the constant and dynamic portfolio allocation and 
hedging problems, where the time-varying covariance matrix was estimated using 
a GJR-AGDCC model, that captures asymmetric volatilities and correlations. For 
the error term we have proposed a flexible infinite mixture of Gaussian distributions, 
which was handled using Bayesian non-parametric approach. We have presented a 
short simulation study that illustrates the differences arising from different assump- 
tions for the errors and shows the adaptability of the DPM model. 

We have employed the proposed approach to solve the portfolio allocation and 
hedging problems using real data of asset returns. In both applications we have 
showed that even though the point estimates for optimal hedge ratios and optimal 
portfolio weights are very similar for Gaussian and infinite mixture models, the 
non-parametric credible intervals are wider. Therefore, the normality assumptions 
forces the investor to be overconfident about her estimates. Moreover, the non- 
parametric model allowed for some one-step-ahead volatilities come from a very 
volatile component, thus making the posterior distribution of covariance matrix 
asymmetric. 

The explained methodology and obtained results are not limited to the two 
specific risk management problems and could be expanded into various other topics 
in applied finance and risk management. 
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Figure 1. True and Estimated Contours of the One-Step- Ahead Returns r t +i 
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Figure 2. Kernel Smoothing Densities of One-Step-Ahead Errors, their Volatilities 
and the Volatilities of the Returns 
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Figure 3. Log-Returns and Histograms 
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Figure 4. The Right and Left Tails of Kernel Smoothing Densities of Predictive 
Returns r t+1 
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Figure 5. The right Tail of the Posterior Distribution of H^ 1 ^ 
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Figure 6. Posterior Distributions of Utility-Based Portfolio Optimal Weights, Ex- 
pected Returns, Variance and Expected Utility 
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Figure 7. Posterior Distributions of GMV Portfolio Optimal Weights, Expected 
Returns and Variance 
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Figure 8. Log-Returns and Histograms of S&P 500 Index and E-mini Futures 
Contracts 
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Figure 9. The Right and Left Tails of Kernel Smoothing Densities of Predictive 
Returns r t+1 
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Figure 10. Kernel Smoothing Densities of the Elements One-Step-Ahead Covari- 
ance Matrix for DPM Model for S&P 500 and E-mini Data 
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Figure 11. Posterior Distributions of Utility-Based Total Hedged Portfolio Optimal 
Ratio, Expected Return, Variance and Expected Utility 
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Figure 13. Kernel Smoothing Densities of Optimal Number of Contracts 
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